Comment on "Nonclassical Smoothening of 
Nanoscale Surface Corrugations" 

In a recent Letter Erlebacher et al. describe the ex- 
perimental observation of nonclassical smoothening of a 
crystalline surface. In their experiment, ripples of typical 
wavelengths Xx — 290 — 550nm in the x-direction were 
formed on Si(OOl) by sputter rippling and then annealed 
at 650 — 750°C. They report that in contrast with the 
classical exponential decay with time, the ripple ampli- 
tude A{t) followed an inverse linear decay of the type 
A{t) — A{0)/[1 + at]. Erlebacher et al. interpret their 
results in terms of a one dimensional step flow model, 
and compare them with theoretical work of Ozdemir and 
ZangwiU (OZ) |]. While we have no reservations about 
the experimental results, we claim that they are incon- 
sistent with the one dimensional step flow model, and 
a fully two-dimensional model is necessary in order to 
describe the behavior of the experimental system. 
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FIG. 1. Amplitude relaxation of an initially sinusoidal grat- 
ing (solid line) measured in units of the height of an atomic 
step, hs. The dashed and dash-dotted lines are fits to an 
exponential decay and an inverse linear decay, respectively. 
Time is in arbitrary units. 

The mass transport rate limiting process for surface 
evolution depends on the typical terrace width, Aa:, and 
on the length £ = Dg/k, where Dg is the adatom sur- 
face diffusion constant and k is the rate of attachment 
of adatoms to step edges. The kinetics of the system is 
diffusion limited (DL) when £ <C Ax and attachment- 
detachment limited (ADL) when £ 3> Ace. Bartelt et al. 
showed 1^ that for Si(OOl) surfaces, £ > bOOnm for tem- 
peratures in the range 670 — 1010°C. This value of £ is 
more than an order of magnitude larger than typical ter- 
race widths in the experiment of Erlebacher et al., and 
therefore their system should exhibit ADL kinetics. 

We studied the relaxation of large-amplitude one di- 
mensional gratings [Q in terms of the same one dimen- 
sional step flow model used by OZ We found that 
unless the system is purely DL, the height of the grat- 
ing decays exponentially in time, agreeing entirely with 
the experiments of Keeffe et al. g . Figure shows that 
even when £ = and the grating amplitude is 50 steps 
(values comparable to the experimental ones), an expo- 


nential decay fits the step model results better than an 
inverse linear decay. Thus, the experimental results of 
Erlebacher et al. are inconsistent with the one dimen- 
sional theory. 

The confusion arises because of the inaccurate state- 
ment made by OZ, claiming to have found inverse linear 
decay in all kinetic limits, but presenting numerical evi- 
dence only for DL kinetics. Our work [Q shows unequivo- 
cally that the decay is exponential in all cases but purely 
DL kinetics (which shows inverse linear relaxation). 

We now argue that a two dimensional model is needed 
in order to describe the experiment of Erlebacher et al. 
The atomic force microscope images in [0 indicate that 
the corrugations are not purely one dimensional. In addi- 
tion to the short wavelength ripples along the a;-direction, 
the surface is corrugated in the transverse y-direction 
with a typical wavelength Xy « 10 A^;. We showed that 
for radially symmetric morphology the diffusion current 
J satisfies J oc dfi/d(^, where fi is the local step chemical 
potential and C — r+2£ h{r)/hs, with hs being the height 
of a step. The generalization of this relation for non ra- 
dially symmetric morphology allows an estimate of the 
order of magnitude of the peak-to- valley directional cur- 
rents Jx and Jy according to Ji cx A/i/ (Ai/2 + 2£Ah/hs) 
for i = X OT y. Here Afi is the peak-to-valley chemical 
potential difference and Ah is the peak-to-valley height 
difference. According to this, Jy/Jx ~ 0.95 in the ex- 
periments of Erlebacher et al. A one dimensional model 
is appropriate only when Jy is negligible compared with 
Jx, which is clearly not the case here. 
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